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ABSTRACT 


This thesis decribes the development of a parallel GaAs 
Meeeeamplifier configuration for the Transient Electro- 
magnetic Scattering Range at the Naval Postgraduate School. 
The implementation of the parallel amplifier configuration 
improved the bandwidth and the signal-to-noise ratio of the 
range. Existing software was used to process signals 
scattered from canonical targets. When interpreted, this 
data was used to optimize the new amplifier configuration. 
The upgraded facility 1s shown to provide measurements 
having excellent agreement with theoretically-predicted 
Signatures of simple canonical shapes. A @emaniftveant 
improvement of the signal-to-noise ratio and bandwidth of 
scattered signatures vis-a-vis that of the previous 
Taboratory configuration is observed. Results’ are 
documented for a target library created to facilitate 


research into a number of scattering problems. 
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A. BACKGROUND 

Electromagnetic analysis and design have tended to be 
historically-based on frequency domain techniques. In this 
domain the impulse response of a target is obtained by using 
a continuous wave (CW) system to sweep frequencies over the 
entire band of interest. Time-domain transient scattering 
measurements offer a faster, less expensive and more 
directly observable alternative. Made possible by the advent 
of fast pulse sources and sampling oscilloscopes, time 
domain techniques use direct pulse excitation to obtain the 
transient scattering response of a target. 

The original Transient Electromagnetic Scattering 
Laboratory (TESL) at the Naval Postgraduate School (NPS) was 
an outside ground-plane range [Ref. 1). This facility, 
initiated in 1980, was limited to symmetric targets due to 
mae requirement for “mirror imaging". In 1983, an indoor 
facility became operational [Ref. 2]. This new facility used 
a shielded anecholic chamber, which eliminated the 
restriction to symmetrical targets as the target was, in 
effect, suspended within the chamber. 

Them LES facilitated research into radar target 
identification using complex pole natural resonances. The 


transient scattering measurements and signal processing were 


originally done using a Tektronix hardware setup. McDaniel 
(Ref. 3] showed good agreement between the measurement of 
impulse responses for a sphere and a thin wire. 

Prior to the current effort, Sompaee [Ref. 4] designed 
and implemented new software for experimental measurement 
control and subsequent signal processing of transient 
Signatures to synthesize scattering impulse responses of 
scale model targets. In addition, the TESL hardware was 
upgraded to improve signal-to-noise ratio. Automatic 
transient scattering measurements were obtained using an 
IEEE bus controller, an IBM PC and a state-of-thewauwe 
digital programmable oscilloscope (DPO). 

Sompaee obtained good agreement between experimental 
measurements and computations for simple canonical targets. 
His noise analysis showed a 10 dB improvement over the 
Tektronix hardware that was used before the upgrade of the 


JWaee be 


5B.  OvVERViny 

The main objective of this research was to incorporate 
two parallel GaAs FET broadband amplifiers to extend the 
frequency range and further improve the signal-to-noise 
ratio (SNR) of the TESL. The hardware layout was optimized 


and the previously existing flexible coax cabling was 


replaced with minimum length hard coax, given the 
restrictions of the existing structure of the chamber. 

Validation of the improved facility was done by again 
comparing experimental results to those theoretically 
predicted for simple canonical targets. A library of high 
fidelity measurements for simple shapes and scale model 
military targets was created. This library is for use in the 
development of new natural resonance extraction algorithms 
and resonance annihilation filters for aspect invariant 
radar target identification. 

This thesis is broken into seven chapters. Chapter II 
describes the TESL and the hardware changes incorporated as 
part of this thesis. The development of the new amplifier 
configuration is summarized. 

Chapter III presents the theory of transient scattering 
measurements. Major time domain facility considerations are 
briefly discussed and the system frequency-domain represent- 
ation, with its mathematical model and transient response 
solution, are presented. 

Chapter IV briefly reviews the acquisition algorithms 
and signal processing developed by Sompaee. It highlights 
how experimental measurements, once processed, were used in 
the development of the amplifier configuration discussed in 


Chapter II. 


Chapter V discusses the experimental results that were 
used to validate the final configuration of the upgraded 
facility. Both sphere and thin wire scattering measurements 
are shown to have excellent agreement with theoretical 
predictions. 

Chapter VI summarizes the target library created. 
Example results from scale model targets are presented. 

Chapter VII summarizes all work and suggests some 


possible continuing research applications. 


Pie ee aiGES [IMPLEMENTED IN THE NPS TESL 


A. LABORATORY DESCRIPTION 

The NPS TESL is located in room 535 of Spanagal Hall. A 
general layout is depicted in Figure 1. The range 
incorporates a metallically shielded anechoic chamber, 6.2 m 
long by 3.1m square. Targets were mounted on a low density 
styrofoam pedestal, 2.18 meters from the frontal plane of 
the antenna for this thesis work. The back wall of the 
chamber iS covered with 46 cm long pyramids while the source 
wall is covered with 21 cm long pyramids. Sides, floor and 
ceiling are covered with longitudinal wedges to direct 
energy towards the absorbing back wall. The interior of this 
chamber, looking towards the source wall, is illustrated in 
Figure 2. Further details of the chamber can be found in 
Pererences 2 and 3. 

Details of the target acquisition hardware and software 
can be found in Reference 4. As part of this effort the 
following hardware was added to the facility: 

- Midisco inductive power splitter; 

~ AVANTEK APT-12066 six-to-twelve GHz amplifier; 
- HP 8349A microwave amplifier; 

iamwedemae) Mil coax cabling; 

- Two cooling fans; 

- One additional horn antenna; and 


- Abacus 25 MHz Cache Memory 80386 Computer with 80387 
and Weitek 3167 Coprocessors 
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Figure 1. NPS TESL general layout 


B. ANTENNA ARRAY 

A 24-inch-square detachable panel was available in the 
center of the source wall of the chamber for mounting the 
antenna array. Two antennas, one transmitting and one 
receiving, were previously used. A new panel for mounting 
three antennaS was needed. Since it was not practical to 
mount all three antennas in a line, the two transmitting 
antennas were mounted above the center line of the panel and 


the receiving antenna below. The three antennas were mounted 





Eanciimne. 2: Interior of anechoic chamber 


to be symmetrical about the center of the panel, as can be 
seen in Figure 3. The antennas were then shimmed to point to 
the position in space where the targets are mounted, at a 
Gistance of 2.18 meters from their face. As will be 
explained in Chapter IV, the metalized mounting panel was 
shielded with the same absorbing material as the source wall 
to eliminate any multi-path effects between the target and 
the mounting panel. Figure 4 shows how the antenna mounting 
panel was masked. 

A shelf was constructed on the exterior of the antenna 
mounting panel in order to mount the amplifiers in close 
proximity to the transmitting antennas. The amplifier 
configuration is depicted in Figure 5. The evolution of this 
configuration is explained in the following sections. Figure 


6 shows the general exterior equipment layout. 


C. AMPLIFIER CHARACTERIZATION 

The first step in the development of the parallel 
amplifier configuration was the characterization of the new 
6-12 GHz GaAs FET unit. The manufacturer’s guaranteed 
specifications for the Avantek APT-12066, one watt 6-12 GHz 
amplifier are listed in Table 1. The device was first 
mounted in a modified heat sink of unknown origin. As the 
device was too wide to fit into the heat sink, the two 


innermost fins were milled out. "The heat Sink was fitteameae 
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Figure 3. Antennas mounted symmetrically 


Figure 4. Masking of antenna mounting panel 
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Figure 6. Chamber exterior equipment layout 
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the length of the device, 7.5 cm. The amplifier was mated to 
the heat sink using type 120 thermal joint compound (EG&E 
Wakefield Engineering). The modified heat sink, with 


amplifier mounted, is illustrated in Figure 7. 


foo 1. MANOUFACTURER’S SPECIFICATIONS FOR NEW AMPLIFIER 


Frequency Noise Power Gain 

Response Gain Gain Figure Output Flatness 

Model (GHz) (dB) (dB) (db) (Watts) (+/-dB) 
Veer Min. Max. Max. Min. Max. 
APT—-12066 One 2s ag, 8.0 ee O Pe 2 


The characterization of the amplifier was done uSing an 
HP 8756A Scalar Network Analyzer and 8350B Sweep Oscillator. 
This was performed in three stages: low input power, medium 
input power and a full specified power test. In all three 
stages, a -20dB HP 33340C attenuator was mounted on the 
output port of the amplifier to protect the HP 11664A 
detector used with the network analyzer. The device was also 
investigated using a bias voltage of less than the specified 
+ 15 volts dc. During testing the device case temperature 
was monitored with a Fluke 80T-150 temperature probe and HP 


3469B multimeter. 
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Figure 7. Amplifier mounting in heat sink 
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At low input power (-30 dBm) it became apparent that the 
the heat sinking described above was inadequate as the 
measured case temperature rose above the specified maximum 
(50 C) after approximately one hour. The device performance 
remained within specification until the case temperature 
rose above 55 C. The bias voltage was decreased to +12 volts 
DC and the case temperature stabilized at 78.4 Cc. At this 
case temperature the average gain across 6-12 GHz was less 
than 26 dB. At a bias voltage of + 15 volts the temperature 
mese above 100 C and did not stabilize. At this point the 
test was terminated for fear of damaging the device. 

Once a small fan waS positioned to blow over the 
amplifier the case temperature stabilized at approximately 
41 C, and this temperature did not vary measurably between 
an input power of -20 or -10 dBm. The device gain remained 
stable at roughly +31 dB across the 6-12 GHz bandwidth. A 
fesponse peak of 33.05 dB was noted at 11.1 +/- 0.05 GHz. 

No measurable difference in the output of the device was 
found when it was biased between +13 and +15 Volts Dc. 
Similar results were obtained when confirming the 
Characteristics of the 1-6 GHz device [Ref. 4]. For all of 
this work both devices were operated with a bias voltage of 
+13 volts DC to keep the case temperatures within the 
mManufacturer’s specification. With the two personal 


computer-type cooling fans used in the final configuration 


LS 


it was found that the device temperature would approach 50 C 
if the room temperature rose above 22 C when biased at +15 
Volts. When biased at +13 C the device temperature would 
stabilize at approximately 44 C with an ambient room 
temperature of 22 C. 

To confirm the device reliability an eight-hour full- 
power test was then conducted. The input power was set at -2 
dBm and the average device response was greater than +31 QB. 
The ambient room temperature at the start of test was + 17.9 
C and rose to + 22.2 C by day’s end. The device response was 
measured periodically throughout the day at 6, 8, 10, and 12 


GHz, aS well as the peak amplification frequency of 11.1 +/- 


of test, the amplifier took approximately one and one-half 
hours to reach a stable temperature. The output of the APT- 


12066 at the end of this test is illustrated in Figure 8. 


D. PRE-AMPLICICATION AND DELAY LINE 

The initial concept for the implementation of the 
parallel amplifiers was to take the output from the pulse 
generator through an inductive power splitter, directly into 
the two GaAs FET amplifiers and on into the transmitting 
antenna, using the minimum possible lengths of hard coax. 
When this was set up on a test bench and fed directly back 


into the sampling head through 40 GB of attenuation, it was 


ie 


TABLE 2. FULL POWER TEST RESULTS 








GAIN 

Case 
Time Temp. 6 GHz 8 GHz 10 GHz 12 GHz 

C dB aB aB aB 
0830 34.0 3 13 9 S09 S S145 30.64 
1005 40.9 ENO SIS: 30:97 31.48 30.40 
ako 5 42.1 30.90 30.96 oe a Os 37 
2 30 40.8 30.96 30.96 Saray, 30.39 
F330 42.8 3.09 1 30.96 31,46 30.36 
1500 43.6 S0787 30.96 31.46 S05 
1625 42.9 30.88 30. 95 Saeed 7 30% 3°60 


found that the the magnitude of the 6-12 GHz pulse was app- 
roximately one-fifth of the 1-6 GHz pulse. Spectral invest- 
igation and simple modeling of the 25 psec rise time step 
pulse from the DPO confirmed this order of magnitude diff- 
erence in the spectral distribution over the 1-6 versus 6-12 
GH2 spectrum. This amplitude difference can be seen in an 
example measurement from the bench testing in Figure 9. 
Figure 9 was plotted directly from the DPoO. The peak 
amplitude of the 6-12 GHz pulse is roughly 12 mVolts. The 


peak amplitute of 1-6 GHz pulse is approximately 60 mVolts. 
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Return signal without pre-amplification 


Figure 9. 


ie, 


This amplitude difference was clearly undesirable, given 
the objective of improving the bandwidth and SNR of the 
TESL. The APT-12066 Amplifier was the highest gain 1 watt 
6-12 GHz amplifier available when purchased. It was decided 
to use an available HP 8349A microwave amplifier to pre- 
amplify the 6-12 GHz pulse. Using the minimum possible 
lengths of hard coax, given the physical chamber structure, 
equipment size and need for air circulation, the pre- 
amplification introduced a 3.95 +/- 0.1 ns delay in the 6- 
12 GHz pulse. This time delay can be seen in Figure 10. 

When this pre-amplifier configuration was implemented, 
with sample measurements taken and processed, (to be 
described in Chapter Three), it was observed that deep 
"valleys" or holes existed in the spectrum in the area of 
spectral overlap between the amplifiers. This observation 
lead to the conclusion of destructive interference between 
the pulses, possibly because of the time delay introduced by 
pre-amplification. It was then decided to delay the 1-6 GHz 
pulse so that the scattered signals arrived at the receiving 
antenna simultaneously. 

To determine the delay necessary a number of scattering 
measurements were made and averaged. The implication of 
these measurement was that a 3.95 +/- 0.1 ns delay was 
required. The +/- 0.1 ns delay was not resolvable by the 


visual images on the monitor or printer being used. 
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Figure 10. 


Inherent noise made it it difficult to measure (or to 
determine on a hard copy) the exact location of the 
subjective front edge of a scattered pulse. 

There was also some uncertainty as to the effect of the 
dielectric fill in the hard coax cable being used) ine 
dielectric constant for polytetrafluoroethylene (Teflon) 
found in handbooks [References 5, 6 & 7) varied from 2.08 to 
2.1. It is suspected that some variance with frequency over 
1-12 GHz exists but no supporting literature was found. 
Using formula 2.1 and the manufacturer-specified inner 
radius ("a"), outer radius ("5") and charactermem ae 


impedance ("Zo") | Re tewwe mee neuse jetta toma dielectric 


constant of 2.0006. [Ref. 9] 
In(b/a) 
(by reek WI JOIN OI (Zeus 


Given the uncertainty in ¢€ it was determined that the 
length of delay line needed was 84 +/- 2 cm, in addition to 
the 11.2 cm feed line used when the measurements were made. 

A number of delay lines were manufactured, varying from 
76 to 102 cm in length. Each was used for a series of 
measurements and the resulting data was processed using a 
deconvolution algorithm and an FFT plotting package [Ref. 
10}. The results were examined. for completeness of their 


spectrum and for evidence of interference in the amplifier 


Ze 


spectral overlap area. It was found that delay lines 
measuring from 92.71 to 99 cm gave good results. The best 
results were obtained using a 97.8 cm delay line. A typical 
scattered signal, obtained using the 97.8 cm delay line 
Semeiguration,; 1S shown in Figure 11. The final delay line 


can be seen in Figure 12. 
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Figure 12. Implementation of 97.8 cm delay line 
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III. THEORY OF TRANSIENT SCATTERING MEASUREMENTS 


A. TIME DOMAIN SCATTERING MEASUREMENTS 

Transient electromagnetic scattering measurements can, 
in essence, be considered an extension of Time Domain 
Refectometry (TDR). In the TESL, the TDR pulse source is 
available from channel 1 of the four-channel test set (HP 
54121A), a sub-assembly of the DPO. This sub-assembly also 
acts as the trigger for the DPO sampling circuits. The rise 
time of the step generator is approximately 25 psec and, for 
this work, waS operated at a pulse repetition rate of 500 
KHz 

The target acquisition algorithm [Ref. 4] is set up so 
that the number of data points in a time series has to be a 
power of 2. Ensemble waveform averaging increases power or 
energy SNR by the number of averages, N. For this work 1024 
data points were acquired in a 20 ns window, averaged 2048 
times. 

The close proximity of the transmitting and receiving 
antennas results in a very large direct coupling of energy 
into the receiving antenna. As the pulse exits the 
transmitting antennas this direct-coupled energy iS many 
orders of magnitude larger than most transient target 
responses. If both were recorded, the desired response would 


be unobservable using the required visual scale factors. 
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This potential problem is avoided by setting the delay 
between the pulse generation and the start of the return 
sampling so as to time window the transient target response. 
With the new amplifier configuration this delay was 40 nsec. 
By that time the residual "ringing" between the antennas was 
maeh smaller than the target response and, as will be 
explained in the following sections, was subtracted out in 


the post-processing. 


lee BeoieoYoTEM REPRESENTATION 

A free-field scattering range system representation is 
shown in Figure 13. The interactions in the scattering 
process within the system are represented by frequency 
domain transfer functions. 

The amplified pulse generator output, X(f), drives the 
transmitting antennas whose combined transfer function is 
femoced by H,~(f). The transmitted field immediately couples 
into the recelving antenna, a) ewe eee nee Urn nicer acts 
Geeen the scatterer, Ho(f), and the absorbing material in the 
chamber, Ho(f). The mutual interaction between the scatterer 
mieeeche Chamber walls is denoted by the two-way arrow, H-(f) 
to Ho(f). The effects from all these sources are summed in 
the receiving antenna, whose transfer function 1s Hp(f). 


Noise is added to the system from thermal emissions within 
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Figure 13. TESL system representation 


the chamber and from the sampling front end of the DPO. 
Total noise is denoted N(f). The net signal sampled by the 
DPO is denoted by Y(f). 

To extract the target Utransfer function, Hao(f),; )ee 


separate measurements are required. The first 1s with no 
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target present and is termed the 
The second is with a target present, 
measurement", 
Lrg 


superscript 


termed the "calibration measurement", 


functions are denoted with a superscript "Cc". 


mee SUmmarizes the transfer functions 


termed the 


used 


"background measurement". 


enGdimatine transrer LunCtCions are denoted with a 
tremchird mrs wren anvecanonltCcade target, 
and the transfer 


The following 


in 


mathematical model developed in the following section. 


TABLE 3. 
Amplified pulse driving the antennas 
Transmitting antennas response 
Receiving antenna response 
Antenna coupling 
Target response 
Chamber response (clutter) 


Interaction between target and chamber 
(including multipath and multiple scattering) 


Target response, target measurement 


Multiple scattering between target and absorber 
target measurement 


Calibration sphere response, calibration 
measurement 


Multiple scattering between calibration sphere 
and absorber, calibration measurement 


Total system noise 


Net signal which is sampled 
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SYSTEM REPRESENTATION TRANSFER FUNCTIONS 


0 i 

Hp (f) 
Hp(f) 
Ha (f) 
Hg (f) 
Ho(f) 


Hso(f) 


Hg (f) 


sea) 
Hg (f) 
HS¢(f) 


N(f) 


Y(f) 


ear Ge e 


C. MATHEMATICAL MODEL 

Each of the three measurements described in Part III B 
can be expressed as a product and sum of their constituent 
transfer functions. | Tmemsesuleann= Signal spectrum at the 


Sampling head for a background measurement will be given by 
¥,(f) = Hp (ft) Hp (fl) Ae) CeCe) eee ee eee (32 


This equation states that the received response to the 
amplified signal is equal to the product of the driver and 
antenna transfer functions times the sum of the antenna 
coupling and chamber response. The system noise is added to 
this product. The following two equations are similar in 
form except that the target or calibration sphere response, 
and their applicable interaction with the chamber, are 


included in the sum. A target measurement is given as 


¥5(f) = Ha(f) Hp(f) X(£) (H,(£) + He(f) + Het) + Heo 


ae) 
(3ieZe 


A calibration measurement is given as 


¥,(f) = Hp(f) Hp(f) X(£) {Ha (ft) + HCCf + Helter 


1 Nee) 
(3s) 


These three equations represent the transforms of the 


measurements obtained directly in the time domain. The post- 


30 


processing described in the next section leads to the 


extraction of the desired transient response. 


D. TRANSIENT RESPONSE POST-PROCESSING 

Deconvolution post-processing makes it possible to 
synthesize the transient scattering response of the measured 
target to any pre-defined incident field whose spectrum fits 
within the pass band of the measurement system. A double 
gaussian pulse was used as the specified excitation in the 
deconvolution process, as will be discussed shortly. With 
the implementation of the parallel amplifiers, the DPO, 
eleritfier configuration and transmitting antennasS are 
roughly matched in the pass band of 1-12.4 GHz. 

To synthesize the target’s transfer function, the first 
Seep 1S to subtract from the calibration and target 
measurements the chamber clutter and the antenna coupling. 
These are embodied in the background measurement, equation 
mere Tie subtraction of equation 3.1 from 3.2 and from 3.3 


leads to equations 3.4 and 3.5. 


¥,(f) = Y5(f) - Y,(f) 


Prgiat. (f) X(t tHe (fs Haa(fy yes Naif) (ae4) 


We(f) = Y(f) - Y,(f) 


II 


Penna (te Ce) (re (yee Bech) lene Ne (f) (355) 


ot 


where N,(f) and Ne(f) represent the difference in the noise 
part of the measurement between the background and the 
target and calibration measurements, respectively. 

The difference of the time domain measurements 
represented by equations 3.4 and 3.5 are then transformed 
into the frequency domain by a fast Fourier transform (FFT). 
An optimal deconvolution estimator, developed by Riad [Ref. 


11), is then formed: 


* 
¥4(f£) Yo(f) 
Xg(£) Hg(£) = --------5-=------ (X9(£) HE (f£)) (3.6) 
Veta eae) se 


The smoothing parameter "C"", establishes a minimum for the 
denominator and prevents noise enhancement at frequencies 
Where the spectrum of Ys approaches zero. The expression 
(Xo (f£) Hg (£) ) represents the computed calibration sphere and 


X5(ft) is the Fourier transform of x5(C), 


X(N) =S{x (t)} (the specified incident pulse) (303 


The double-gaussian pulse used in this work is specified 


bys 


Z 2 
-act—-t ) -a t=} 
x (V=c @ a oe e al o) (3 39 
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Cy and C, are chosen such that 
and 


+ 00 
[x at = 0 (3.10) 


— 69 


The double gaussian waveform of equation 3.8 will have two 


specified 10% pulse widths, At,, such that 





(8) | 
ere Oot lore Rat 2 ia) 

It can be shown that the estimator of equation 3.6 
provides the best least-squares fidelity for the deconvolvcd 
target response. Equation 3.6 represents the frequency 
domain scattering response for the specified incident pulse 
whose Fourier transform is X)(f)- Two programs written by 
Morgan are used to numerically generate the theoretical 
response, (X9(f) eae) Nie of the calibration sphere [Ref. 2”. 
Called Mie and SPRSCT, the latter uses magnitude and phase 
data generated by a Mie series in the former to generate the 
transient response of the sphere due to the double gaussian 
incident pulse. 

The subtraction and optimal deconvolution routines are 
part of the software developed by Sompaee [Ref. 4]. In the 
next chapter it will be shown how the results of these 
algorithms were used in the development of the amplificr 
configuration. 
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iV.|)POsT-PROGESe ate 


A. DECONVOLUTION ALGORITHM 

The subtraction and optimal deconvolution algorithm, as 
described mathematically in the previous chapter, are 
embedded in the TESL computer software support system. This 
program waS originally written by Morgan for the first TESL 
Tektronics microcomputer [{Ref. 2]. It was converted to 
Fortran 77 by Sompaee and installed in the microcomputer 
used in the TESL [Ref. 4]. The program has graphic 
Capabilities and allows the operator to obtain a hard copy 
at each step of the deconvolution procedure. The data files 
used in all the TESL software support systems have a 
standardized format. 

The program starts by asking the operator for the name 
of the file containing the theoretically computed transient 
scattering by the calibration sphere target. This file is 
produced by successive runs of two programs MIE.BAS and 
SPRSCT.FOR. After reading the file, the program displays 
the time-series plot of the theoretical scattering. The 
calibration sphere scattering is then transformed into the 
frequency domain and displayed so that the bandwidth can be 
checked. The bandwidth of the theoretical scattering 
obtained from the Mie series and SPRSCT programs must not 


exceed that of the measured data. Frequency and time-domain 
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plots for the 15.1 cm calibration sphere used in this work 
are shown in Figure 14. 

The next step of the program is to read in the files for 
the measured calibration sphere and the measured background 
data. Both files are overlaid on the same axis and displayed 
Mmemeene operator. At this point the operator has the 
capability of shifting one of the waveforms to the left or 
right if he observes that time shifting has occurred. This 
capability was not found to be needed during this work. An 
example target and background measurement comparison is 
shown in Figure 15. The background waveform is then 
subtracted from the calibration measurement and the 
difference is transformed into the frequency domain. 

The program then asks the operator for the target file 
name and its background file name. The operator has the 
option of using the background file subtracted from the 
calibration measurement or a different background file. In 
each case, over the course of this work, a background file 
was used which was obtained either immediately before or 
after the target measurement being processed. The program 
then follows the same procedure described for the 
Salibration sphere: overlay, option to shift in time, 
subtraction and transformation. An example subtracted 


waveform is shown in Figure 16. 
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Figure 14. Time and frequency domain plots, theoretical 
calibration sphere 
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Figure 16. Subtracted waveform 
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The three frequency-domain transforms, the FFTs of the 
numerically computed calibration sphere, the measured 
calibration shpere and the measured target response are then 
formed into Riad’s optimal deconvelution estimator given in 
equation 3.6. The program then asks the operator for the 
smoothing parameter "C", This parameter is kept as small as 
possible to prevent oversmoothing and loss of spectral 
content. It was found that C = 0.009 gave good results and 
was used throughout this work. The frequency domain 
estimator formed by Raid’s method yields ae best 
approximation for the scattered signal produced by the same 
incident plane wave used for the computed calibration 
sphere. This frequency domain estimator is displayed to the 
operator. The program then transforms it into the time 
Gomain and the result is again displayed. The operator has 
the option of saving the deconvolution result, reprocessing 
with a different smoothing parameter, or processing another 


target. 


B. USE IN AMPLIFIER-ANTENNA DESIGN AND VALIDATION 

In the example subtracted waveform, as presented in 
Figure 16, the residual system noise both before and after 
the scattered signal can be observed. This noise iS random 
and has no apparent harmonic content. When trial 
measurements were first made with the new amplifier-antenna 


configuration, a seemingly harmonic content in the noise was 


So 


observed. This harmonic content was observed only in a 
target-subtracted waveform. It was not observed in the 
background waveform or when a background measurement was 
subtracted from another background waveform. It was 
concluded that the harmonic content could be from residual 
multi-path scattering between the metalized antenna mounting 
panel and the target. When the panel was masked, as can be 
seen in Figure 4, the harmonic content of the noise in the 
subtracted waveform disappeared. 

The frequency domain estimator formed by the post- 
processing provided the first clue to another potential 
problem: possible destructive interference between the two 
pulses after the 6-12 GHz pulse had been delayed in time by 
pre-amplification. In Figure 17, a deep hole can be seen in 
the frequency domain estimator at roughly 5 GHz. When the 
raw data from the same measurement was transformed into the 
frequency domain, using an FFT, several nulls were observed 
in the area of spectral overlap between the two amplifiers. 
A number of trial measurements were made and the same 
observation was made in the FFT spectrum and in the 
frequency domain estimator between 5 and 6.5 GHz. 

To confirm what was actually being transmitted into the 
chamber, each antenna lead was in turn disconnected and 
terminated through two 20 dB attenuators and a 50 ohm 


terminator. A series of scattering measurements were made 
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with each antenna terminated. These measurements were then 
processed using the deconvolution algorithm followed by FFT 
transformation. In all cases the spectral content from 5 to 
6.5 GHZ was found to begrien. 

The deconvolution algorithm was also used to determine 
the length of delay line required to bring the two pulses 
into phase with each other at the receiving antenna. A 
number of trial measurements using different-sized spheres 
were made with each length of delay line. The measurement 
data was then put through the post-processing and the 
resulting frequency domain spectral estimator was compared 
to the theoretically predicted spectral plot. Each delay 
line was also lengthened by approximately 1 cm through the 
addition of a male-male and female-female 141 mil coaxial 
Section. The measurements and subsequent processing were 
then repeated. Thus, by trial and error, the final delay 
line length was determined. 

In the next chapter the deconvolution results for three 
different canonical targets will be presented. These were 
obtained using the final amplifier configuration. In 
addition, an SNR calculation based upon measurements will be 


Considered 
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V. SYSTEM VALIDATION 


A. CANONICAL TARGETS 

To validate the fidelity of the final equipment 
configuration the time domain results of the deconvolution 
program are compared with numerically generated waveforms 


for simple canonical targets. Three targets were used: 


-~ A 12.28 centimeter diameter metal sphere, 
- A 8.10 centimeter diameter metal sphere, 


- A 10 centimeter long, 2.36 millimeter diameter thin wire. 


Three programs, MIE, SPRSCT and TDIE-DG, written by 
Morgan [Ref. 2], were used to generate the theoretical 
waveforms for the targets. MIE computes the frequency domain 
Magnitude and phase for the metallic spheres. These are then 
used as inputs to SPRSCT to compute the transient scattering 
response to the specified double-gaussian incident waveform. 

TDIE-DG was used for the thin wire calculations. This 
program solves a time-domain integral equation to generate 
the induced current on the wire for a specified double- 
gaussian pulse incident plane wave field. The scattered 
field is then calculated by a numerical space-time 
integration using the induced currents. The input parameters 


for these programs are summarized in Table 4. 
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TABLE 4. NUMERICAL COMPUTATION 


Sphere radius in meters, 12 cm sphere 
Sphere radius in meters, 8 cm sphere 
Bistatic angle in degrees 

Scattering plane 

Time window in ns 

Narrow 10% pulse width in ns 

Wide 10% pulse width in ns 

Wire length in meters 

Wire radius in meters 


Scattering angle in degrees 


Scattering angle in degrees 
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INPUTS 


0.061423 


0.04050 


0.00118 


90 
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The numerically generated theoretical waveform and the 
measured deconvolution waveform for the two spheres and the 
thin wire are overlaid in Figures 18, 19, 20 and 21. Thin 
wire results are presented for both 90 and 45 degree 
incidence angles. The time base has been shortened as 
compared to Figure 14 so that the detail can be examined. It 
can be seen that the agreement between the plots is 
excellent. The high level of fidelity of the TESL and the 
accuracy of the deconvolution algorithm are apparent from 
these overlays. 

Some residual noise and clutter can be observed as a 
Slight ripple on the measured waveforms. These signal 
pollutants appear to be less than what was seen in previous 
work. The system noise estimation presented in Reference 4 
was repeated for comparative purposes and is presented in 


mae next section. 


Eee oYoLEM NOISE ESTIMATION 

In order to compare the system signal-to-noise ratio to 
that calculated during previous work on the TESL [(Refs. 3 & 
4), the same thick wire target was again measured. This 
target is a 10 cm long, 5.2 mm diameter, copper wire. The 
measured response of this wire 1S shown in Figure 23. Two 
consecutive background measurements were also made for the 
purpose of the noise estimation. One background measurement 


was first subtracted from the target waveform. The 
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subtracted signal waveform was then time-windowed to include 
only the signal duration. This time series is denoted by 
x (Cs) ine mewo consecutive background measurements were then 
subtracted from each other to eliminate chamber clutter. 
This record is denoted as N{(t_). Statistical jeunmarte aie 
these two records are found in Tables 5 and 6. The signal- 
to-noise ratio (SNR) was then calculated as 


THE 80] Sd) 


Silos : : 

ne = 

I N XN (t,) 
Nn 








(525) 


The resulting SNR obtained uSing equation 5.1 was 26.6 
Glist When this effort started, it had been expected to 
replace the one watt 1-6 GHz amplifier with a two watt 
device. The two watt amplifier is unavailable until 1990 
and a further increase of SNR is expected when the new 
amplifier is installed. The thick wire measurement was 
repeated several times and similar results were obtained. 
Sompaee (Ref. 4] reported a SNR of 27.31 dB. When his target 
measurements were re-examined, it was found that the 
seemingly harmonic content of the noise existed, which lead 
to the masking of the antenna mounting panel in this work. 
Thus, his SNR may have been contaminated with a multi-path 
contribution which was eliminated during this work. As well, 


it was found that by allowing the amplifiers to reach 


on 0) 


thermal stability before measurements were taken produced 
better results. For this work, the hardware was allowed to 
warm up for approximately two hours before measurements were 
taken. The SNR calculation was also repeated uSing the 
target waveforms discussed in the next chapter. Direct 
comparisons could then be made using the same target, same 
number of ensemble averages, and the same software used for 
the computations. These calculations demonstrated the 
improvement achieved in the SNR. Some target aspects 
demonstrated a very large response at frequencies greater 
than 7 GHz. These records showed an improvement in the SNR 
of up to 8 dB. In no cases examined was there a decrease in 
the SNR as compared to the target records from the single 


amplifier configuration. 


Slt 








eel Ss sotATISTICS OF THICK WIRE SIGNAL PLUS NOISE 


WAVEFORM 
MEAN S-4132525 
S’D DEVIATION 5.03768 =4 
MEDIAN 3.6864E-5 
MAXIMUM 1465 = 3 
MINIMUM =e JOO => 


TABLE 6. STATISTICS OF NOISE WAVEFORM 


MEAN lL -409660E509 
oa. DEVIATION 25 524E=05 
MEDIAN 958108 S06 
MAXIMUM tee 2OS Li OS 
MINIMUM Sry ae ESD 1 Os 
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VI. TARGET LIBRARY 


The fidelity of the parallel amplifier configuration was 
demonstrated uSing the: simple canonical targets discussed in 
Chapter V. The scattering response of the canonical targets 
is numerically computable to a high level of accuracy. This 
is not the case for more complex scattering targets. A 
primary purpose of the NPS TESL is to support research into 
aspect-invariant radar target identification, based on 
natural resonances. To facilitate this research a library of 
measurements on scale model aircraft targets was assembled. 
Six 1/72 scale-model targets of military aircraft were 
selected. Each plastic model has been carefully coated with 
Silver paint to excite a scaled transient scattering 
response, duplicating the aircraft modeled. 

Each target file contains 14 measurements. Each 
measurement is constructed from 2048 ensembled averages of 
1024 data points in a 20 ns window. The same measurements 
were made for each of the aircraft targets as follows: 


- A calibration sphere measurement 


A background measurement 
- Target measurement, 30 degrees from nose-on (30°) 
- A second background measurement 


- Target measurement, nose-on (0°) 


54 


-~ A third background measurement 
- Target measurement, broadside (90°) 
- A fourth background measurement 


~ Target measurement, tail-on (180°) 


A fifth background measurement 


- A target measurement, viewed from the top with the 
wings horizontal 


A sixth background measurement 


- A target measurement viewed from the top with the 
wings vertical 


~ A seventh background measurement. 


It took roughly one hour to acquire each target waveform and 
its background pair. The DPO acquired the data in nine 
Minutes per measurement, after which it took roughly 15 
Minutes for the IEEE 488 bus transfer to the computer. Hard 
copies of each target response were then made before 
proceeding to the next measurement. Figures 24 to 27 are 
examples of the waveforms obtained. Each has been time 
windowed from the original 20 ns record so detail can be 
seen. Table 7 lists the full-size dimensions of the 


aircraft modeled. 


oye) 


TABLE 72 £ULi 


Target number 
Overall length 
(meters) 


Overall height 
(meters) 


Wingspan (meters) 


Tailplane span 
(meters) 


SIZE DIMENSIONS OF TARGETS RECORDED 


12.20 16.00 19.43 16.94 153503) 2072 


By iS us) 4.80 Sj ey) 4.51 5s OS 4.46 


10.96 13.295 13.05 11.43 10200) eae 


Unknown. $5275 Sie: 649.2 5258 N/A 
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VII CONCLUSIONS 


A. SUMMARY 

The thrust of the effort described in this thesis has 
been the development, implementation and validation of a 
mae lel amplifier configuration for the NPS TESL. The 
resultant increase in bandwidth and SNR will enhance the 
guality of the measurements taken to support a number of 
research applications. To support those efforts, an 
extensive library of high fidelity target measurements for 
both canonical and complex scattering targets was created. 

The first stage of this effort was the development of 
the parallel amplifier configuration. The TESL software 
Support system was used extensively in this development. A 
scalar network analyzer was used to characterize the new 6- 
12 GHz amplifier and the first parallel configuration tested 
led to frequency domain modeling of the relative magnitudes 
of the spectral content of the pulse generator from 1-12 
GHz. This modeling lead to the preamplification of the 6-12 
GHz amplifier which, in turn, lead to the delay of the 1-6 
GHz pulse. The delay was required to eliminate destructive 
interference between the pulses caused by the time-delay 
introduced by the preamplification. 

The parallel amplifier configuration was then validated 


by COmparing the outputs of the deconvolution post- 
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processing for canonical targets to the numerically computed 
scattered waveforms due to a double-gaussian incident pulse. 
Very good agreement was shown, which attests to the accuracy 
of the experimental procedure and the deconvolution 
algorithm. A noise analysis was presented and the causes of 
the SNR improvement seen in the new hardware configuration 
were examined. 

The final stage of this work was the creation of a high 
fidelity library of scattering responses for canonical and 
complex targets. The library has been described and example 


waveforms presented. 


lee FUTURE CONST DERATIONS 

In Chapter V, Section B, the prospect of furan 
increasing the SNR by replacing the one watt 1-6 GHz 
amplifier was mentioned. Other possible hardware 
improvements include the use of a matching two watt 6-12 GHz 
amplifier, should one become available. In addition, the 
preamplifier currently employed could be replaced with a 
more compact low-noise device. This latter measure would 
facilitate rack mounting the time domain equipment, DPO, DC 
power supply, sampling front-end and amplifiers. This 
equipment could then be more readily removed and replaced 
with alternate hardware configurations for use in obtaining 


broadband scattering measurements. 


6zZ 


Under consideration is the use of a more mature and more 
widely-used technique employing a  stepped-frequency 
continuous wave system. It would use a network analyzer 
front-end as a coherent receiver. Frequency domain 
measurements could then be made and directly compared to the 


library compiled in the time domain. 
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